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Aims. To investigate the very high energy (VHE: >100 GeV) 7-ray emission from the high-frequency peaked BLLac 1ES 0229+200. 
Methods. Observations of 1ES 0229+200 at energies above 580 GeV were performed with the High Energy Stereoscopic System 
(HESS) in 2005 and 2006. 

Results. 1ES 0229+200 is discovered by HESS to be an emitter of VHE photons. A signal is detected at the 6.6a level in the HESS 
^ . observations (41.8 h live time). The integral flux above 580 GeV is (9.4+ 1.5 s tat ± 1.9 sys t) x 10 -13 cm _2 s _1 , corresponding to 
~1.8% of the flux observed from the Crab Nebula. The data show no evidence for significant variability on any time scale. The 
observed spectrum is characterized by a hard power law (F = 2.50 ± 0.19 s tat ± 0.10 aya t) from 500 GeV to ~15 TeV. 
Conclusions. The high-energy range and hardness of the observed spectrum, coupled with the object's relatively large redshift 
(z = 0.1396), enable the strongest constraints so far on the density of the Extragalactic Background Light (EBL) in the mid- 
infrared band. Assuming that the emitted spectrum is not harder than Fi nt ~ 1.5, the HESS data support an EBL spectrum oc A -1 
and density close to the lower limit from source counts measured by Spitzer, confirming the previous indications from the HEGRA 
data of 1ES 1426+428 (z = 0.129). Irrespective of the EBL models used, the intrinsic spectrum of 1ES 0229+200 is hard, thus 
locating the high-energy peak of its spectral energy distribution above a few TeV. 

Key words. Galaxies: active - BL Lacertae objects: Individual: 1ES 0229+200 - Gamma rays: observations 



1. Introduction 

The active galactic nucleus (AGN) 1ES 0229+200 was 
initially discovered in the Einstein IPC Slew Survey 
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(Elvis et al. 1992) and later identified as a BL Lac ob- 
iect (|Schachter et al.l 1993). It is now classified as a high- 
frequency peaketil3LJLac_(HBL) due to its X-ray-to-radio 
flux ratio (Giommi et al. 1995) , The HBL is hosted by an 
elliptical galaxy with absolu te magnitude Mr — —24.53 
(Falomo & Kotilai nen 1999)[ located at a redshift of z = 
0.1396 (Woo et al. 2005)[ Based on its spectral energy dis- 
tribution (SED) IBS 0229+200 is suggested as a potential 
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source of VHE 7-rays ( Stecker, de Jager fc Salamon 1996 
Costamantc & Ghisellini 2002). However, despite several 
attempts, it has not been previously detected in the 
VHE regime. The Whipple (Ide la Calle Perez et al. 2003| 



IHoran et al. 20041). HEG RA |( Aharonian et al. 2004a)[ and 
Milagro [(Williams 2005) collaborations have each reported 



upper limits on the flux from 1ES 0229+200 during various 
epochs. The most constraining upper limit (99.9% confi 
dence level) on the flux is I(>410 GeV) < 2.76 x 10 
cm -2 s _1 , based on ~1 hour of HESS observations in 2004 
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(Aharonian et al. 2005) The present discovery of VHE 7- 



rays from 1ES 0229+200 makes it the most distant ob- 
ject detected at multi-TeV energies. As a result its VHE 
spectrum can provide important information on the EBL 
(Hauscr & Dwck 2001) in the mid-infrared band, where few 
measurements exist. 



2. Observations and Analysis Technique 

1ES 0229+200 was observed with the HESS array of i mag- 



ing atmospheric-Cherenkov telescopes (Hinton 2004) for a 
total of 70.2 h (161 runs of ~28 min each) in 2005 and 
2006. After applying the standard HESS data-quality se- 
lection, 98 runs remain yielding an exposure of 41.8 h 
live time at a mean zenith angle Z m eap = 46°. The data 



are calibrated as detailed in Aharonian et al. (2004b) and 
the standard HESS analysis tools (Benbow 2005) are used. 
The eve nt-selection criter ia are performed using the stan- 
dard cuts (Benbow 2005) resulting in a post-analysis energy 
threshold of 580 GeV at Z mcan . A circular region of radius 
cut = 0.11° centered on 1ES 0229+200 is used for the on- 
source data. The background (off-s ource data) is estima ted 
using the Reflected- Region method (Berge et al. 2007) As 
the source direction was positioned ±0.5° relative to the 
center of the field-of-view of the camera during the ob- 
servations, this method allows the simultaneous estimation 
of the background using the same data as the on-source 
measurement. The significance of any excess is calculated 



following the method of Equation (17) in Li & Ma (1983) 
The energy of each event passing selection is corrected for 
the absolute optical efficiency of the system, using effi- 
ciencies determined from simulated and observed muons 



(Aharonian et al. 2006b) This corrects any potential long- 
term variations in the observed spectrum and flux due 
to the changing optical throughput of the HESS system. 
Results consistent with those presented are also found us- 
ing independent calibration and analysis chains. 



3. Results 

The results of the HESS observations of 1ES 0229+200 in 
2005 and 2006 are shown in Table [TJ A significant excess 
of 261 events (6.6<r) from the direction of 1ES 0229+200 is 
detected in the total data set. The on-source and normal- 
ized off-source distributions of the square of the angular 
difference between the reconstructed shower position and 
the source position (9 2 ) are plotted in Figure Q] for all ob- 
servations. The background is approximately flat in 9 2 as 
expected, and there is a clear excess at small values of 9 2 
corresponding to the observed signal. A two-dimensional fit 
of the observed excess finds the shape to be characteristic of 
a point source, located at (aj2ooo = 2 h 32 m 53.2 s ± 3.1g tat ± 
1.3| yst , feooo = 20°16'21" ± 44£ at ± 20^). The excess, 




0.15 
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Fig. 1. The distribution of 9 2 for on-source events (points) 
and normalized off-source events (shaded) from observa- 
tions of 1ES 0229+200. The dashed line represents the cut 
on 9 2 applied to the data. 
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Fig. 2. The VHE spectrum observed from 1ES 0229+200. 
The line represents the best \ 2 fit of a power law to the 
observed data. Only the statistical errors are shown. 



named HESS J0232+202, is consistent with the position 
(jSchachter et all 1993) of the blazar (a J20 oo = 2 h 32 m 48.6 s , 
<5j20oo = 20°17'17") as expected, and is therefore assumed 
to be associated with 1ES 0229+200. 

Figure shows the time-average photon spectrum ob- 
served from 1ES 0229+200. The best X 2 nt of a power 



law (dN/dE ~ E r ) to these data yields a photon index 
T = 2.50 ± 0.19 sta t ± 0.10 syst , and a x 2 of 3.6 for 6 degrees 
of freedom. No evidence for significant features, such as a 
cutoff or break, are found in the measured spectrum. 

Assuming the determined photon index of V = 2.50, 
the observed integral flux above 580 GeV is I(>580 GeV) 
= (9.4±1.5 st at±l-9 S yst) x 10~ 13 cm~ 2 s -1 . This corresponds 
to ~L8% of I(>580 GeV) measured by HESS from the 



Crab Nebula (Aharonian et al. 2006b) Table Q] shows the 
flux measured during each year of HESS observations. All 
values are below the previously published VHE upper limits 
mentioned earlier. The \ 2 an d degrees of freedom for fits 
of a constant to the integral flux versus time when binned 
by months within each year and years within the total, are 
also shown in Table [1] As the \ 2 probability for each fit, 
as well as for fits of the flux binned nightly within each 
dark period, is greater than 0.2, there is no evidence for 
variability on any time scale within the HESS data. 
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Table 1. The MJD of the first and last night of HESS observations of 1ES 0229+200, the live time of the observations, 
the number of on- and off-source events measured, the on/off normalization (a), the excess, and the significance of the 
excess are given. In addition, the integral flux above 580 GeV (assuming T — 2.50), and the corresponding percentage of 
the Crab Nebula flux above 580 GeV are shown. The x 2 , degrees of freedom (NDF), and % 2 probability, P(x 2 ), for a fit 
of a constant to the flux binned by dark period within each year, or yearly within the total, are also given. 



Epoch 


MJD 

First 


MJD 
Last 


Time 

w 


On 


Off 


Q 


Excess 


Sig 

M 


I(>580 GeV) 
[10- 13 cm^s- 1 ] 


Crab 

% 


X 2 , NDF 


P(x 2 ) 


2005 


53614 


53649 


6.8 


246 


2238 


0.09160 


41 


2.7 


6.8±3.1 s tat±1.4 syst 


1.3 


1.6,1 


0.21 


2006 


53967 


54088 


35.0 


1344 


12304 


0.09136 


220 


6.1 


10.0+1. 7 st at±2.0 syst 


1.9 


1.5,3 


0.68 


Total 


53614 


54088 


41.8 


1590 


14542 


0.09140 


261 


6.6 


9.4±1.5 sta t±1.9 sys t 


1.8 


0.8,1 


0.37 



4. Discussion 



VHE photons from extragalactic sources are expected 
to suffer absorption from interactions (7 VHE 7 EBL — > e + e~; 
Gould & Schroder 1967| with EBL photons (see, e.g., the 
review of Aharonian 2001). The energy-dependent opacity 



E (TeV) 



t(z, E) causes a deformation of the observed VHE spec- 
trum (F b s (E) = Fint(-E') e _T ^ ,£; ') which depends on the 
EBL SED (Figure [3]). With assumptions about a source's 
intrinsic VHE spectrum, limits on the EBL density can be 
derived. Recently, the hard 0.2—2 TeV spectra observed 
from 1ES 1101-232 and H 2356-309 constrained the EBL 
around 1— 2/im to a significantly low level, close to the 
lower limits given by the integrated light from resolved 
galaxies { A haronian et al. 2006a, Aharonia n et al. 2007[) . 



Apart from 1ES 1426+428 (| Aharonian et al. 20030 . 
1ES 0229+200 is currently the only VHE source at z>0.1 
whose spectrum has been measured up to 10 TeV. 
Since the peak of the 7-7 cross-section occurs at A* w 
1.4 {E 1 / 1 TeV) /im, these VHE spectra probe the EBL in 
the near (NIR) to mid- infrared (MIR) range (« 2— 20 /mi), 
where direct estimates are missing due to increasing 
foreground emission from interplanetary dust and where 
EBL models differ significantly in the slope of the SED 
(Figure [3]). This difference is exemplified by comparison of 
the models by Primack et al. (2005, hereafter PrimOS) and 
Stecker et al. (2006, hereafter SteckOd). Aside from nor- 
malization, the NIR-MIR sp ectrum in the Prim05 model, 
as well as in the models 



by Kneiske et al. (2004) 
the characteristic decline of the old stellar component 
in the SED of galaxies, which behaves approximately as 
vl(v) oc A -1 at these wavelengths. The corresponding pho- 
ton number density, n(e) cx e , makes the optical depth 
almost constant with energy between one and several TeV 
( Aha ronian "20011 . This causes a characteristic flattening 
feature in the attenuation curves. This flattening will be im- 
printed in the observed VHE spectra for an intrinsic power- 
law spectrum. In contrast, the EBL decline is much flatter 
than A -1 for both the baseline and fast evolution versions 
of the Steck06 model. This yields a rapidly increasing opac- 
ity above 1 TeV which should translate into a strong and 
continuous steepening of the observed VHE spectra (unless 
counter-balanced by upturns in the intrinsic spectrum). It 
should be noted that, regardless of the EBL model, the op- 
tical depth is large (i.e. r > 1) for all energies relevant here. 
Therefore, the observed spectrum of 1ES 0229+200 is very 
sensitive to t(E) and hence to the shape of the EBL SED. 




A[yu.m] 

Fig. 3. SED of the EBL at wavelengths most relevant for 
the VHE spectrum of 1ES 0229+200. The upper axis shows 
the energy scale corresponding to the peak of the 7-7 cross- 
section. Data compilation from Aharonian et al. 2006al (see 
references therein). Open symbols: Lower limits given by 
the integrated light from resolved sources. From left to 
right: Hubble data (jMadau fc Pozzetti 2000)), Spitzer data 
(jFazio et al. 2 004), ISO data (tri angles; fAltieri et al. 19991 
lElbaz et al.~20 02). Spitzer data (|Dole et al. 2006|) . Dashed 
lin es: Fast evolution (u pper) and baseline (lower) models 
of iStecker et al. (2006)1 Dot-dashed line: Model [PrimOS) 
of |Primack et al. (2005)] Solid lines: Scaled (xl.6) Prim05, 
modified above 2 jam by a range of allowed MIR slopes, 
i.e. those which yield an intrinsic spectrum rj nt = 1.5 and 



follows 1.15 » 1.5-0.23. 



0.1 S y S t. The fluxes at 10 /im are 2.15 



and 3.1 nW m 2 sr 1 , respectively. 



The intrinsic spectrum of 1ES 0229+200 can be recon- 
structecQ by multiplying the individual data points by 
e r( z ,E) For the twQ steck06 modelfl the 1ES 0229+200 
absorption-corrected spectra are well fit by power-law func- 
tions (dN/dE ~ E~ Fint ) with hard photon indices Tint = 
0.6 ± 0.3 (baseline EBL) and Tint = 0.1 ± 0.3 (fast evolu- 
tion EBL). Such hard spectra are not easily explained by 
commonly used leptonic or hadronic scenarios for the 7- 



1 Throughout the paper, a cosmology with Ho=70 km/s/Mpc, 
n m =0.3 and Q,a=0.7 is used for the calculation of t(z,E). 

2 The optical depths were calculated by integrating over the 
respective EBL SED. The methodol ogy was checked with the 



prescription 
redshift (z— 
the intrinsic spectrum 



in [Stecker et al. (2006) for the nearest tabulated 
0.117), obtaining consistent values (Ar <0.05) for 
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ray emission of BL Lac objects, unless invoking extreme as- 
sumptions like mono-energetic particle distributions, either 
in the shocked plasma (see, e.g., |Katarzynski et a~2 006) 
or as a cold wind with very high (—10') Lorentz factors 
(| Aharonian et al.~2 002). Further, a specific dependence of 
source parameters on redshift would be required, in order 
to explain the absence of these hard features in all nearby 
(z<0.1), less-absorbed VHE-bright HBL. 

The Prim05 model yields r int = 1.92 ± 0.22. Since this 
model falls below the recent lower limits given by Spitzcr 
source counts at MIR ( Fazi o et al. 20041 see Figure[3]) wave- 
lengths, and thus underestimates the attenuation at several 
TeV, modified EBL shapes consistent with the lower limits 
from source counts are considered in the following. 

In standard blazar models, VHE 7-rays are produced 
by particles accelerated in shocks. For a wide range of 
conditions the res ulting photon spectrum is expected to 
have Tint > 1.5 (| Aharonian et al. 2006a[) . This can be 
used to constrain the EBL. It should be noted that 



Stecker et al. (2007) have demonstrated that intrinsic spec- 
tra as hard as Fj nt ss 1.0 are possible in relativistic shock 
acceleration scenarios. In such situations, the derived EBL 
SED constraints will be relaxed. Following Aharonian et 
al. (2006a), the intrinsic photon spectrum is reconstructed 
from an EBL SED shapqfl and then tested for compatibil- 
ity with r^t > 1.5. Different EBL shapes can be derived 
for examination, either by scaling a fixed EBL shape or 
by changing the NIR/MIR flux ratio (thus testing different 
EBL slopes). 

In the observed energy range of 1ES 0229+200, a sim- 
ple scaling of EBL SED shapes « A -1 , like the P1.0 curve 
in Aharonian et al. 2006a or the Prim05 shape, does not 
strongly affect the reconstructed intrinsic spectrum due to 
the weak energy dependence of the optical depth. Indeed, 
the spectrum is only modified due to the difference in r 
for low and high energy 7-rays, which is not the case here 
where the EBL mainly causes attenuation. The limit for 
the overall scaling factor from the HESS 1ES 0229+200 
data is 0.8x P1.0 (i.e. P0.8), or 1.8x Prim05. The limit on 
the flux is above those obtained from other VHE sources 
(| Aharonian et al. 2006a|) . 

Different EBL spectral slopes are tested using an EBL 
template constructed from a fixed scaling of the Prim05 
shape up to 2/im and an interpolation above this wave- 
length. The MIR-EBL slope is varied using a linear (in log- 
log scale) interpolation between 2, 10, and 30/im, where 
only the 10 /xm flux is changed (see Figure [3]). The scaled 
EBL level is chosen to be 1.6 x Prim05, which corresponds 
to the EBL upper limit at 1—2 /jm as derived by Aharonian 
et al. (2006a). This NIR EBL level allows for the highest 
possible MIR fluxes, since a high flux at 1— 2/im softens 
the intrinsic spectrum of 1ES 0229+200 in the energy band 
measured here by HESS. Thus, much higher overall EBL 
levels are compatible with Tint > 1.5. For a flux at 10 /im 
of 2.15 nW m -2 sr _1 the intrinsic spectrum is Tint = 1.5. 
A limiting 10/zm flux of 3.1 nW m~ 2 sr _1 is derived by 
taking into account the errors on the photon index (i.e. 



3 In the calculation of r for these test-SEDs, the effects of 
galaxy evolution were approximated by modifying the cosmo- 
logical dependence of the photon number density as (1 + z) 3 ~ k . 
With k=1.2, the values of t(z) are in good agreement with those 
from the Prim05 model at z < 1. Regardless, evolution effects 
are small at z — 0.140. The difference between including or not 
this effect is Ar < 0.08. 
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Fig. 4. The broad-band SED of 1ES 0229+200. The ATOM 
optical measurements (stars) are simultaneous to some of 
the 2006 HESS data. All other data are non-simultaneous. 
Both the observed and intrinsic (i.e. de-absorbed with a 
1.6 x Pmm05 EBL shape, yielding r int = 1.58 ± 0.22) 
VHE spectra are shown, each with the two highest- 
energy HESS points combined for visibility purposes. 
Using a lower EBL density results in a softer intrin- 
sic spectrum, but even the very low density of Prim05 
yields Tint < 2.0. The X-ray data are from Einstein 
(2 keV), ROSAT (1 keV; IBrinkmann et al. 1995)) , and 
BeppoSAX (butterfly; IDonato et al. 2005)) observations. 
The radio, NIR (2MASS) and optical data are from 



the NED archive, Urry et al. (2000) 



Perlman et al. (1996) 



and lSchachter et al.K 1993). The ATOM fluxes and the data 
marked with open triangles do not have the contribution of 
the host galaxy subtracted. 



where r int = 1.15 « 1.5 - 0.23 stat - 0.1 syst ). The EBL 
SED between 2 and 10 /xm (oc A~") is thus constrained to 
a slope a > 1.1 ± 0.25. A T^t ~1.5 spectrum is also ob- 
tained with an EBL SED at the level of the source counts, 
connecting the optical data points with the Spitzer data 
at 8//m (for an a ~ 0.9). Clearly the HESS data do not 
support flatter EBL shapes in the MIR, confirming a trend 
already suggested by t he HEGRA data on 1ES 1426+428 
(] Aharonian et al. 2003)) . 



5. Conclusions 

The HESS spectrum of 1ES 0229+200 leads to two impor- 
tant conclusions. First, unless the intrinsic BLLac spec- 
trum is considerably hardei@ than 1.5, the EBL density 
in the MIR range follows a spectrum cx A -1 and is very 
likely close to the lower limits given by galaxy counts as 
determined by Spitzer and ISO data. This implies that the 
sources resolved by Spitzer seem to account for the bulk of 
the diffuse 5— 10/im background. Second, almost indepen- 
dently of the details of the EBL SED, the intrinsic spectrum 
of 1ES 0229+200 is significantly hard up to -10 TeV. Thus 
the high-energy peak of its SED (shown in Figure 2]) is lo- 
cated at energies above a few TeV, among the highest ever 
seen from a blazar. The extension of the SED to VHE en- 
ergies provides clear evidence that non-thermal processes 



4 See, e.g., Stecker et al. (2007) for scenarios that demonstrate 
Tint as hard as 1.0 is possible in shock acceleration models. 
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are at work in 1ES0229+200. The SED also includes opti- 
cal (Bessel B and R filters) measurements of 1ES 0229+200, 
rriR = 16.2 and tub = 17.2, conducted with the ATOM tele- 
located at the HESS site. These 
in some of the 2006 HESS obser- 



scope (Hauser et al. 2004) 
data are simultaneous wit 



vations. The ATOM fluxes arc constant in time and fall 
below archival values. Although only limited evidence for 
a VHE-optical flux correlation exists |(Albert et al. 2006) | 
the low optical state might suggest a correspondingly low 
VHE state. These results strongly motivate further VHE 
observations of 1ES 0229+200, as well as contemporaneous 
observations at lower energies to enable SED modeling. 
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